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The quasi-two-dimensional modeling of the small adiabatic per- 
turbation on the background of the stationary configuration of the 
1/^ I selfgravitating gas with the weak transverse nonhomogeneity approx- 

^^ ■ imation is presented. The space periodic character of the solution of 

^ . this system is proofed. 

1j : 1 The Physical Model 



An effects arising during an expansion of the plane stationary adiabatic wave 
with a small amplitude on the background of the stationary configuration of 
the self gravitating gas have a volume nonlinear character. These effects 
r> I arising in the process of competition of the two basic mechanisms: the self 

c^ I gravitating and thermal expansion of particles. Three physical parameters 

define the solution of this problem: $o - the characteristic gravitational po- 
tential, To - the characteristic temperature for the thermal expansion and 7 
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- the index of adiabatic. The stationary configuration of the self gravitating 
gas has cyhndrical symmetry [|l|. We use the term "radius" to fix the bound- 
ary between the areas of the gas and vacuum. The wave under consideration 
perturbed to a small extent the initial distribution of the speed, density and 
radius. Smallness of the perturbation is of principle character. The mathe- 
matical description of this model is represented by the system of equation of 
motion and continuity (three dimensional nonstationary partial differential 
equations), Poisson equation and algebraic equation of state: 

dv VP 

^ + (zJ-V)t; = -^-V$; (1) 

|^ + V-(pt;) = 0; (2) 

V^$ = 4nGp; (3) 

P = Ap', (4) 

there P is pressure, v - thermal speed, $ - gravitational potential, p - gas 
density, t - time, A = const. This system describes the ideal classic gas with 
self gravitation. In the system of coordinate related to the stationary wave 
the derivative of time disappear [|I| from the system (^). An addition to (^) 
we'll use also the radius component of the equation of motion multiplyed by 
p{z, r)r: 

1 dvl dP d^ 

2 The Relations of the Quasi-Two-Dimentional 
Averaging Method 

The QTDA method is applied to finite three dimensional space problem 
with separate direction. A rigorous and consistent quasi-two-dimensional 
approach must based on the cross-section averaging of three dimensional 
physical model equations. This procedure allows one to proceed from three- 
dimensional distributions of physical quantities (such as potential, mass den- 
sity etc.) to their one- dimensional averaged distributions [Q. The following 



relations of the QTDA method have been used: 

Siz) L(z) 

S{z) L{z) 

I V. i(r,.)d. ^ ^^{{Mz)ySU)} + / i^d. (8) 

where n - external (as referred to the section S{z) ) normal to the side surface; 
■^ is the angle between the axis r and the normal n to the side surface of the 
cylinder; 
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dz 
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S{z) = / V<l>{r,z)-ds; (9) 



S(z) 



{^{z)fS{z) = I <^{r,z)ds, (10) 

S{z) 

where index ^ marks the cross section averaging values. All characteristics 
of the three dimensional problem are described in the new equations by the 
relation establishment between new functions - the average quantity variable 
of physical value and the boundary meanings of the initial function. In this 
formulas the boundary meanings are represented by the contour integrals. 
It is seen that all non-one dimensional effects are located in the contour 
integrals. 

3 The QTD Model of the Weak Transverse 
Nonhomogeneity Approximation 

The QTDA equations are the result of application of the QTDA method 
represented by (3) - ([T0|) to equations of continuity, Poisson, and (2) in the 
system of coordinate related to the wave. In this model we studing the effects 
related to the small perturbations of the physical values. According to this 



the smallness of following parameters have been supposed: 

d ^ (^'f-yr , (12) 

Po 

where cq and po are the sound speed and density on the axis z in the nonper- 
turbated part of the cylinder. The smallness of these parameters correspond 
to the weak transverse nonhomogeneity approximation. We study the equa- 
tions of continuity, Poisson, state and (2) in the 0- approximation of small 
parameters (|T^). The structure of the QTDA equations depended from con- 
tour integrals, which are simplified in the 0-approximation because in this 
case: 1) cos^ = 1; 2) n is parallel to V$; n ■ Cz = 0. Therefore only the 
Poisson equation has a nontrivial summand after the averaging. Thus the 
QTDA equations under consideration become as follows: 

kvr = -^ipr-'-^'+c, (13) 

z 7 — i 

where C = const; equation of continuity: p^v^ = const; Poisson equation: 

^^!^-A$-(.)=4.Gp-(.). (14) 

It is not possible to obtain the second summand from one dimensional model. 
The averaging of the 0-approximation of the equation (2) represent the de- 
pendence of the radius of the cylinder from the coordinate z: 

- = MR) = v.K(a)-, (15) 

where a is the surface density, K{a) = 2(1- ) '^'^''^~''^ ■ Let's estimate benefit 
of transverse motion energy in the full energy of the system. For this purpose 
we'll substitute v = Vr + Vz iia ( [T3| ) and compare the energy of longitudinal 
and transversements we find out, that the latters are unessential far from the 
wave front Az, 

Vz{R)K{a) 



Az ^ max 



-Aa 



(16) 



where F^ - one of the set from ({v^^f,^{p^y'\2 ^^ ,2|C|j. Aa - is 

a changing of surface density along the length Az. It means that volume 
effects are more valuables than valuable surface ones. We'll use in following 
the averaging values without index ^. 

4 The QTD Description of the Self Gravitat- 
ing Gas 

For future investigation the dimensionless variables (the coordinate, density, 
pressure and potential) have been introduce: 

^j PO -MD •i'o 

where Lj - is the Djeence's length, cq - sound speed. There are characteristic 
values marked with index "0" in this formulas. The dimensionless equations 
are: 

^ = u{l--f5^+^), (19) 

u=iip + 6)il--fS^+^), (20) 

6 = (3uP, (21) 

where ' = ^, d^ = {^ — 'j)S'^~^^dT. The meaning of density related to sound 

speed is (5s =( - J , so the (^1]) is not equal to infinity. The first integral 
of (gl|) is 

ip = gid)+C, (22) 

where C = const. From the system of equations 

^ = (l-^5^+i)(^(5) + 5 + C), (24) 

^ = /?.^^ (26) 



the second integral of the system ( pS] ) has been received and has a form of 
"energy conservation \aw":E = ^ + U, where u plays a role of speed, E is 
the constant "energy", and the "potential" 

U = -l[C + g{5)r + ^{-^ + 5'^). (26) 

The local minimum of the potential U correspond to finite solutions of the 
system (E^). It is according to result have received in |1| for 7 = 2. 
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